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We demonstrate Mach-Zehnder-type interferometry in a superconducting flux 
qubit. The qubit is a tunable artificial atom, whose ground and excited states 
exhibit an avoided crossing. Strongly driving the qubit with harmonic exci- 
tation sweeps it through the avoided crossing two times per period. As the 
induced Landau- Zener transitions act as coherent beamsplitters, the accumu- 
lated phase between transitions, which varies with microwave amplitude, re- 
sults in quantum interference fringes for n — 1 ... 20 photon transitions. The 
generalization of optical Mach-Zehnder interferometry, performed in qubit 
phase space, provides an alternative means to manipulate and characterize 
the qubit in the strongly-driven regime. 
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The development of artificial atoms with lithographically defined superconducting circuits 
presents a new paradigm of quantum solid state physics [1J, allowing the realization and ex- 
ploration of new macroscopic quantum phenomena j2l|3l|4l|5l|6l[ZllEl|9), an d a ls° holding 
promise for applications in quantum computing ifTOll . Of the various effects demonstrated with 
qubits, of greatest importance are time-dependent coherent phenomena. Those include the ob- 
servation of Rabi oscillations in charge, flux, and phase qubits 0161 El El ED, entanglement 
of two ifTTTl and three lfT2l qubits, coherent oscillation lfT3l and bifurcation ifLU in multilevel 
systems, and the demonstration of basic elements of coherent control lfT31 fT6l fTTll . Artificial 
atoms strongly coupled to photons have opened the arena of "circuit quantum electrodynamics" 



Here we demonstrate an application of superconducting qubits to quantum physics, realized 
in a strongly driven flux qubit, and described in terms of a Mach-Zehnder (MZ) interferometer. 
The conventional MZ setup uses two beamsplitters: the first divides an optical signal into two 
coherent waves which travel along paths with different effective lengths, and the second re- 
combines and superposes these waves, leading to quantum interference fringes in the measured 
output signal. In a driven qubit, following an idea discussed in Ref. [20|, the beamsplitters can 
be realized by Landau-Zener (LZ) transitions at a level avoided crossing. Over one oscillation 
period of the driving field, the qubit is swept through the avoided crossing twice (Fig. 
Starting from the marker, at the first LZ transition (time ti), the ground state |0) is split into 
a coherent superposition of the ground and excited states, |1) and |0), which, after evolving 
independently and accumulating a relative phase A6'i 2 , interfere at the second LZ transition 
(time t 2 ). The corresponding qubit-state energy evolution (first period, Fig. [TJB) between the 
recurrent LZ transitions (shaded region) provides a phase-space analog to the two arms and the 
beamsplitters of an optical MZ interferometer (top left, Fig,[TjB). The interference phase 
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depends on the magnitude of the qubit energy detuning excursion for times t\ < t < t%. The 
interference fringes in the occupation probability correspond to integer and half-integer values 
of A9 12 /2n. Known as Stiickelberg oscillations with Rydberg atoms ET1I221 . this mechanism 
has application to quantum control ll23l . 

The qubit MZ intereferometer differs in a number of ways from an optical interferometer. 
Firstly, instead of a photon, the interferometry is performed using the quantum state of a qubit. 
Secondly, in the qubit, we have the interference of paths in phase space rather than in coordinate 
space; the phase A#i 2 (EqC]) is determined by the qubit level splitting, which plays the role of 
the optical path length. Finally, being more fragile than photons and easier to decohere, qubit 
states can be manipulated in a coherent fashion only at relatively short time scales. 

We use a periodic driving signal, a harmonic variation of the qubit detuning e(t): 

H = -- (Aa x + e(t)a z ) , e(t) = e + A rf cosut, (2) 

with A the tunnel splitting, e the detuning proportional to dc flux bias, and the rf field 
amplitude proportional to the rf flux bias (see supplementary material). In this case (Fig[l)B) 
we have cascaded LZ transitions which occur when the driving amplitude exceeds detuning, 
giving rise to the inteference fringes at > \s \ (Fig|l£). While the phase A6>i 2 equals 
the shaded area in Fig. [TjB and is dependent on A Y f, the total phase gained over one period, 
9 = j^§e{t)dt = 2ireo/hu, equals the difference of the shaded and unshaded areas and is 
independent of A x ^. As consecutive pairs of LZ transitions (consecutive MZ interferometers) 
interfere constructively when 6 = 2nn, the fringes will appear around the resonance detuning 
values 

£o,n = nhv, n = 0,l,2,..., u = uj/2n. (3) 

Another interpretation of this condition is that the sequential LZ transitions excite multiphoton 
resonances. 
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While coherent multiphoton resonances between discrete states of an rf-driven charge qubit 
have been reported ll24l Q, and multiphoton transitions used to drive Rabi oscillations in a 
flux qubit |25j|26l, in these works, as well as in the earlier work on quantum dot systems ll27l 
1281 . only a few-photon transitions could be observed, with coherence quickly weakening at 
increasing rf amplitude ll29ll . In contrast, we were able to observe coherent resonances of very 
high order, up to n = 20, which requires driving the system at a high rf amplitude. The fringes 
for high n are as clear as those for n ~ 1, indicating that the qubit preserves a substantial 
amount of coherence even in the strongly-driven regime. 

We realize a tunable artificial atom with a niobium persistent-current qubit (Fig. a 
superconducting loop interrupted by three Josephson junctions POl . When the qubit loop is 
threaded with a magnetic flux fq « $o/2, the system exhibits a double- well potential-energy 
landscape (Fig. |6j). The classical states of the wells are persistent currents Iq with opposing 
circulation, described by energy bands ±e /2 = ±/q$ 5/q linear in the flux detuning 5fq = 
fq — 3*0 /2- The double- well barrier allows quantum tunneling of strength A, opening the 
avoided crossing near 5f q = (Fig. UK)- Detuning the flux tilts the double well and, thereby, 
modifies its eigenenergies and eigenstates. The qubit states are read out using a DC-SQUID, 
a sensitive magnetometer that distinguishes the flux generated by the circulating currents. The 
device was fabricated at MIT Lincoln Laboratory (see supplementary material). 

We drove transitions between the qubit states by applying a 1 [is rf pulse (Fig. |2jB) at fre- 
quency v and rf-source voltage Vrf OT1 . After a short (~10 ns) delay, we read out the qubit state 
by driving the DC SQUID with a 20 ns "sample" current 1$ followed by a 20 fis "hold" current. 
The SQUID will switch to its normal state voltage V$ if h > -^ sw o (-^s > -^ S w 1-*' corres P on ding 
to qubit state |0) and |1). By sweeping the sample current and flux detuning while monitoring 
the presence of a SQUID voltage over many trials, a cumulative switching-distribution function 
is generated, revealing the "qubit step" (Fig. |2|C). At specific values of flux detuning, the rf 
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field at v = 1.2 GHz becomes resonant with the energy level separation, allowing n-photon 
absorption, Eqj3l this results in a partial population transfer between the qubits states, manifest 
as regularly spaced "spikes" in Fig. [2jZ. We obtain ID scans of the "switching probability" P sw 
(the population of state |0)) shown in Fig. [2)D by following a flux-dependent sample current 
J sw o < Is < I sw \ (dash-dotted line in Fig. Et). Such ID scans are then accumulated as a 
function of the rf source parameters (Fig . O and v (Fig.©. 

The switching probability Psw (color scale) vs. qubit flux detuning 5fq and voltage 
at frequency v = 1.2 GHz is shown in Fig. |3] (see supplementary material). The n-photon 
resonances, labelled by n = 1 ... 20, exhibit MZ interference fringes (I. . .VI) as a function of 
Vrf. As we show below, the fringes exhibit a Bessel-function dependence, J n (A), so we call the 
step-like pattern in Fig.|3]a "Bessel Staircase". For each of the n-photon resonances, we took 
a higher-resolution scan (e.g., Fig. |3] inset) and fitted the resonance areas and widths in Fig. 
(see supplementary material). 

Multi-photon transitions at the resonances (Eq.|3l) in the strong driving regime, \A r f\, hu ^> 
A, occur via fast LZ transitions. The notion of quasistationary qubit levels ±~(£q + A 2 ) x l 2 is in- 
adequate in this regime and, instead, we use a different approach, transforming the Hamiltonian 
(Eq.0! to a nonuniformly rotating frame, H = e -(Wy.W' e (W(*)'« , where <j>(t) = A sin cut 
with dimensionless rf field amplitude 



The rf field disappears from the detuning term, reappearing as a phase factor of the off-diagonal 



nhv ~ e we can replace the phase factor e l ^ by its n-th Fourier harmonic, J n (X)e %nu 
with J n the Bessel function. The resulting effective Hamiltonian is of a "rotating field" form: 



(4) 



term: H 



Eqo z + Ae l ^a + + h.c). Given that A ^ hu, near the n-th resonance 




(5) 
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with A„ = AJ n (X). The resonance approximation (Eq. O) describes transitions at an arbitrary 
ratio Arf/hv. Standard Rabi dynamics analysis of the Hamiltonian © with the initial state 



|0) gives the time-averaged occupation probability of the excited state = ^\A n \ 2 /((sq — 
nhu) 2 + | A n | 2 ). This expression predicts Lorentz- shaped resonances of width 5e = \ A n \. The 
result, a sum of independent contributions with different n, 

1 |A n | 2 

2 „ (£o - nhvy + |A n | 2 

is displayed in Fig. [if. The agreement with the observed resonances is striking: the oscillations 
in rf power, described by J n (A), accurately predict both the overall profile of the fringes (Fig.© 
and the fine details, such as positions of the nodes (Fig©. 

In the frequency dependence of P S w for voltages V r f = 71 mV r ms an d Vrf = 7.1 mV rms 
(Fig. |4j), the resonances approach the qubit step as frequency decreases, in accordance with the 
linear energy vs. flux-detuning dependence. MZ interference fringes are again visible. The 
number of resonances increases at low frequencies, due primarily to the frequency dependence 
of A and, in lesser part, a frequency-dependent mutual coupling. 

Our analysis of peak profile accounts for the relaxation and dephasing, as well as for the in- 
homogeneous broadening due to low-frequency noise. These effects can be separated from one 
another by considering the peak areas A n which, in contrast with the widths of the resonances 
w n , are not affected by inhomogeneous broadening. The standard Bloch approach yields ll32l 



7^ , ufc = -L + + 1 (7) 



with T12 the longitudinal and transverse relaxation times, and T 2 * describing the inhomogeneous 
broadening. These are aggregate relaxation times averaged over the periodic qubit detuning, 
which, in the operating limit e ^> A, tends to overestimate Ti and underestimate T 2 compared 
with their values at the avoided crossing. 



Fig. |5K shows the Bessel dependence of the n = 1 and n = 5 resonance areas fit by 
Eq. |7] including (blue) and omitting (red) times Ti 2 . The corresponding resonance widths and 
their fittings are shown in Fig. |5(B. Fig. |5t and Fig. [5)D show the resonance area and width 
respectively for 10 resonances, including n — 19. Fitting the areas and widths yields self- 
consistent estimates: Ti » 20 fjs,T 2 » 15-25 ns,T 2 * « 5-10 ns, and A/h « (2tt)4MHz. The 
Ti, T 2 , and A estimates are similar to those reported in Ref. [26|. The nearly linear behavior at 
the nodes of J n [see Fig.|5K] indicates that the decoherence is small compared with the splitting: 
TiT 2 (A//i) 2 fa 250 for n = 1, and decreases slightly for n = 5. The fit/data discrepancy for the 
first fringe for n — 1, which disappears as n increases, is traced to ~ 20% thermal population 
of the excited state due to its proximity to the qubit step (see supplementary material). 

This MZ interferometry technique finds application to qubit characterization and model val- 
idation, two increasingly important research areas in quantum information science. In addition 
to coherence times, which can be obtained by multiple means, MZ interferometry allows the 
direct calibration of the microwave amplitude driving the qubit through the Bessel argument A; 
we found the rf mutual coupling to be Mq = 100 ± 2 fH over all 20 resonances. The agree- 
ment between the two-level Hamiltonian in Eq. Eland the observed resonances n = 1 ... 10 
in Fig. |3] is remarkable. The MZ technique also reveals shortcomings of the two-level model 
at strong driving. For example, the influence of a second MZ interferometer at the avoided 
crossing between the first and second excited states results in the moire-like pattern observed 
for resonances n > 12. We also note here an observed (~ 0.1 GHz) shift in the resonance 
positions at strong driving [Fig. |3] inset]. Both effects require the presence of higher-excited 
states modelled by the full qubit Hamiltonian ll30l l26ll . The high stability and coherence of 
the strongly driven qubit, even at n = 20 photon transitions, illustrates not only the potential 
for nonadiabatic control methods [23|, but also indicates the high potential of niobium devices 
fabricated in a fully-planarized, scalable process for superconductive quantum computation. 
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Supplementary Material 



Device fabrication and parameters The device was fabricated at MIT Lincoln Laboratory on 
150 mm wafers in a fully-planarized niobium trilayer process with critical current density Jq ~ 
160 A/cm 2 . The qubit's characteristic Josephson and charging energies are £j ~ (2nh)300 
GHz and Eq rs (27rfr)0.65 GHz respectively, the ratio of the qubit JJ areas is a ~ 0.84, and 
the tunnel coupling A fa (27r/i)0.005 GHz. The qubit loop area is 16 x 16 /mi 2 , and its self 
inductance is Lq « 30 pH. The SQUID JJs each have critical current I c q ps 2 /xA. The SQUID 
loop area is 20 x 20 /im 2 , and its self inductance is £§Q ~ 30 pH. The mutual coupling 
between the qubit and the SQUID is M ~ 25 pH. An experimental spectroscopic analysis of 
the first five energy levels of this qubit can be found in Ref. fl26l . 

Implementation The experiments were performed in a dilution refrigerator with a 12 mK 
base temperature. The device was magnetically shielded with 4 Cryoperm-10 cylinders and 
a superconducting enclosure. All electrical leads were attenuated and/or filtered to minimize 
noise. 

The data were obtained at a repetition rate of 200 Hz, allowing approximately 5 ms for qubit 
equilibration (initialization) between trials. In Fig|3l 3000 steps were taken along the flux axis 
and 250 steps were taken along the amplitude axis. Each point comprised an average of 100 
trials. The data set was accumulated continuously over approximately 4.5 days. The Fig.|5]data 
set was obtained from 2D scans acquired at the same step size in flux and amplitude, but over 
extended amplitude ranges for higher-n transitions and with each point averaging 2000 trials. 
Although the onset of the 20-photon resonance is difficult to observe in Fig.|3]due to the color 
shading, we clearly observed it when scanning to higher amplitudes. The plotting and fitting 
algorithms used to obtain Fig. |3] and Fig. |5] were fully scripted and identically applied to each 
resonance. 
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Temperature effects The thermal energy in frequency units is k^T/h ps 1 GHz as deter- 
mined from the qubit step in Fig. 2D. The qubit is dc -biased, prepared, and measured far from 
the avoided crossing, e ,n ~ hv n > k^T ^> A. The qubit can be operated in this limit, because 
the equilibrium population is predominantly in the ground state (particularly for n > 3), and 
the coherence time scales are much longer than the rf drive oscillation period over which the 
Mach-Zehnder physics occurs. Where residual excited state population does exist, e.g., n = 1 
due to its proximity to the qubit step, it acts to offset the Bessel function dependence of the 
spectroscopy peaks. Since population inversion should not occur for our qubit in the limit of 
steady-state spectroscopy, the first lobe height is clamped at ~ 50% (in our experiment, the 
other lobes do not exceed 50% even with the offset). When the offset is removed and the Bessel 
fit is made, the data/fit discrepancy is observable for the first lobe of n = 1. The effect is less 
pronounced for n = 2 and n = 3 as they are further from the qubit step with a smaller and 
smaller equilibrium excited-state population. 

Differentiation from previous work The multiphoton transitions observed here are distinct 
from Tien-Gordon photon-assisted transport 1 33 1 or above-threshold ionization of atoms (34|, 
both of which utilize a continuum of final states. Furthermore, the Bessel dependence is a 
quantum-coherent interferometric effect and, thereby, differs from the classical, non-linear 
phase-locking mechanism that may also lead to a Bessel dependence, e.g., as with Shapiro 
steps [35 1 . The Landau-Zener and Mach-Zehnder physics observed here is a coherent interfer- 
ometric effect and, thus, differs from previous Landau-Zener work in the incoherent limit ll36ll . 

In our previous work in Ref. [ 26], we observed single-photon, two-photon, and three-photon 
transitions by weakly driving the qubit at frequencies u — 1 ... 20 GHz. In the present work, we 
are able to observe up to n = 20 photon transitions by strongly driving the qubit at relatively low 
frequencies v = 0.2 ... 1.4 GHz, allowing sufficient "headroom" for multi-photon transitions 
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to high-order between the ground and first-excited states. 
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Figure 1 : Mach-Zehnder interference in a strongly driven qubit. (A) Starting at the dot marker, 
the qubit state is swept by an rf field. After a LZ transition at the first avoided crossing (time ti), 
the resulting superposition state of |0) and |1) (dashed lines) accumulates a phase A#i 2 (shaded 
region), and interferes at the return LZ transition (time £2)- The qubit state is subsequently 
driven away from the avoided crossing and then returns to the starting flux position. This 
single period of qubit evolution is a single MZ interferometer. Depending on the interference 
phase A0i2, amplitude may build in the excited state. (B) The corresponding qubit energy 
variation induced by a periodic rf field, Eq|2l results in an equivalent optical cascade of MZ 
interferometers (MZ#l-#3, top) with resonance condition Eq. |3] (C) The population of the 
qubit excited state, Eq0 as a function of rf amplitude and detuning £ - Note the interference 
fringes (I, II, III, ...) at > e and the multiphoton resonances at e = nhv. 

Figure 2: Multiple resonances in a strongly-driven flux qubit. (A) Circuit schematic of the 
three-junction flux qubit (inner loop) with circulating current Iq and the DC SQUID readout 
(outer loop); Josephson junctions are indicated with an x. A time-dependent flux /(£) oc e(t) 
threading the qubit is a sum of the flux bias due to the dc current I b and a pulsed ac current 
at frequency v irradiating the qubit and driving transitions between its quantum states. The 
SQUID is shunted by two 1 pF capacitors to lower its resonance frequency. Resistors mark the 
environmental impedance isolating the SQUID. (B) The time sequence for the rf pulse (duration 
1 \x% and rf-source voltage V T f) and SQUID sample current J s . A repetition period of 5 ms allows 
for equilibration between trials. (C) A cumulative switching-probability distribution of the qubit 
as a function of Is and the qubit flux detuning 5fq under rf excitation at ^0.12 Vmis an d 
v = 1.2 GHz. Multi-photon transitions are observed between the qubit states |0) and |1) and 
are symmetric about the qubit step (5fq = m$ ). (D) The ID switching probability P sw 
extracted from Fig.Et (dash-dotted line scan). 

Figure 3: Multi-photon interference fringes: a "Bessel- staircase". Switching probability Psw 
is plotted as a function of qubit flux detuning 5fq and voltage at frequency v — 1.2 GHz. 
n-photon resonances are labelled 1 ... 20. Each n-photon resonance exhibits oscillations in 
P sw due to a MZ-type quantum interference that results in a Bessel dependence J n (\) with A 
the rf amplitude scaled by hu, Eq.@. Roman numerals mark the interference fringes of J n (X) 
(solid white lines). The n-photon resonances are symmetric about the qubit step (0 m$ ). Inset: 
close-up of the n = 4 photon resonance. 

Figure 4: Frequency dependence of the multiphoton interference fringes. Qubit switching prob- 
ability Psw is plotted as a function of frequency v and flux detuning 5fq in the limit of (A) 
strong driving (V^ = 71 mV rms ) and (B) weak driving (V^ = 7.1 mVnns)- The n-photon 
resonances approach the qubit step as frequency decreases. Symmetric patterns in the J n (A) 
peaks and valleys due to the MZ-type quantum interference are clearly observed. 
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Figure 5: Analysis of the resonance area and width. (A) Resonance area A n versus voltage 
for the n = 1 and n = 5 photon transitions. The Bessel-dependence J n (\) is observed over 
several lobes. The data are best fit by functions that include decoherence (blue line) rather than 
omit it (red line). As illustrated by the insets, decoherence becomes more pronounced as photon 
number increases. (B) The resonance width w n versus voltage for the n = 1 and n = 5 also 
exhibits a Bessel dependence. The area (C) and the width (D) plotted for resonances n — 1 . . . 9 
and n = 19. 



Figure 6: Qubit double-well potential and the corresponding energy band diagram for the two- 
lowest qubit states. The circulating currents of the qubit correspond to the states |0) and |1) of 
the uncoupled left and right wells. These single- well states are coupled with tunneling strength 
A by the potential barrier; the eigenstates of the coupled system are denoted "e" and "g". 
Varying the flux detuning 5f q tilts the double-well potential and, thereby, modifies the energy 
band structure of the qubit. Far from 5f q = 0, the eigenstates and eigenenergies resemble the 
single-well states and energies. Near 5f g = 0, the system exhibits an avoided crossing. 
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